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ABSTRACT 

We have searched three hot molecular cores for submillimeter emission from the nucleic 
acid building-block pyrimidine. We obtain upper limits to the total pyrimidine (beam- 
averaged) column densities towards Sgr B2(N), Orion KL and W51 el/e2 of 1.7 x 
lO^^cm"^, 2.4 X lO^^cm"^ and 3.4 x 10^"'cm~^, respectively. The associated upper 
limits to the pyrimidine fractional abundances lie in the range (0.3 — 3) x 10~^°. 
Implications of this result for interstellar organic chemistry, and for the prospects of 
detecting nitrogen heterocycles in general, are briefly discussed. 

Key words: astrobiology - ISM: individual (Orion KL, Sgr B2(N), W51 el/e2) - 
ISM: molecules - line: identification. 
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1 INTRODUCTION 

Molecular clouds contain many organic molecules that 
are known to be important in biochemistry. Astronom- 
ical observations, particularly at radio wavelengths, al- 
low us to determine the chemical composit ion and char- 
acter i stics of this molecular inventory (e.g . iDickens et alJ 
1200 ll: ICharnlev. Ehrenfreund fc KuanI lioOlfl . Following in- 
corporation into protostellar disks and comets, these mole- 
cules, or their descendants, were probably the major 
source of volatile organic material available to the early 
Earth iChvba et alj iTQQCf) . Studies of molecular cloud 
composition therefore enable us to quantitatively address 
the issue of the connection between interstellar chem- 
istry, the organic composition of primitive Solar System 
material, and the origin, evolution and distribut i on of 
Life in the Galaxy (e .g. lEhrenfreund fc CharnlevI l200Ct 
lEhrenfreund et alj|20^ . 

Many organics that are known, or strongly suspected, 
to be present in interstellar clouds, are fundamental compo- 
nents of the large organic macromolecules that are central 
to biochemistry. Examples of these are sugars and amino 
acids, the respective building blocks of polysaccharides and 
proteins. The simplest of these, glycolaldehyde and glycine, 
are both identified in the interstellar medium jHoUis et alJ 
l2000l:lKuan et alJl2003all . 

Until recently, definitive detections of interstellar 
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ring compounds have been scarce (e.g. ethylene ox- 
ide, iDickens et alJ IIQQTT) and previous searches for im- 
idazole, cyanoform, pyrrole and pyrimidine were unsuc- 
cessfulJSimoii fc SimonI Il973l: iMvers. Thaddeus fc Link^ 
119801: llrvine et alJIlQSlIT Recently we have tentatively de- 
tected the azaheterocyclic compounds 2H-azirin e (C-C2H3N) 
and aziridine (C-C2H5N) at mm wavelengths llKuan et alJ 
|k)03b; Charnlcv, Ehrenfreund & Kuan 20^. Two (differ- 
ent) tentative lines of Aziridine have also been claimed by 
IDickens et al.l i200l|) . These observations suggest that bio- 
chemically important ring molecules may await detection. 
Of these potential discoveries, a key interstellar molecule for 
Astrobiology, comparable in importance to glycine, would be 
pyrimidine (C-C4H4N2), the unsubstituted ring analogue for 
three of the DNA and RNA bases: thymine, cytosine and 
uracil. Interstellar pyrimidine was unsuccessfully se arched 
for at 46 GHz 30 years ago JSimon fc SimonI Il973^ . Evi- 
dence for various purines and pyrimidines in space, includ- 
ing pyrimidine, comes from the fact th at they have been 
detected in meteoritic organic matter jStoks fc Schwart^ 
1981, 198^, and also may be components of C omet Hal- 
ley's CHON dust jKrueeer. Korth fc Kisseilll99J) . 

The molecular composition of hot molecular cores 
is known to largely refiect the solid state chemistry 
that occurred on grain surfaces, prior to their depo- 
sition into the gas after pro tostellar dust heating (see 
lEhrenfreund fc CharnlevI l2000h . Reactions between HCN 
and VyCN (CH2CHCN) on the surface of dust grains has 
been s uggested as a possible source of interstellar pyrim- 
idine dSimon fc SimonI 119731) . Observations of hot cores. 
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Table 1. Source List 



Source 


R.A.(B1950) 


Dec(B1950) 


Vlsr 








(km s"-'-) 


Orion KL 


05'' 32'" 47^00 


-05° 24'30".0 


8.0 


Sgr B2(N) 


Ijh 44m lo^20 


-28° 2l'l5".0 


64.0 


W51 el/e2 


ig*" 21™ 26^30 


+14° 24' 39". 


60.0 



known to contain other rings believed to have formed on 
dust prior to evaporation, such as C-C2H4O ( ethylene oxide) 
Dickens etaPllQQTl: iNu mmeli n et"all Il998l) and C-C2H3N 
Charnley. Ehrenfreund & KuanI boM) . do show the very 
high abundances of the HCN and VyC N molecules re- 
quired for formation of N-heterocycles ( Ziur ys &: Turner! 
1 19861: iTurnedlTggil : iNummelin fc Bergmanlll999l) . Alterna- 
tively, experiments indicate that energetic processing of ices 
(e.g. AUamandola, Bernstein & Sandford 1997) leads to an 
organic residue c ontaining many polycyclic aromatic hydro- 
carbons (PAHs) IGreenberg et alj.2000..) . One may speculate 
that similar processing of interstellar ice analogues contain- 
ing ammonia, molecular nitrogen, and hydrogen cyanide, 
may similarly produce N-heterocycles like pyrimidine. 

Hence, searches of hot molecular cores, employing high- 
quality spectroscopic data, may offer the best chance of a 
pyrimidine detection. In this paper we report the results 
of a search towards three massive star-forming regions for 
submillimeter emission from pyrimidine. 



2 OBSERVATIONS 

We carried out a program of submillimeter observations us- 
ing the James Clerk Maxwell Telescope (JCMT)^ on Mauna 
Kea, Hawaii during Semester OIB over the period September 
9-13, 2001. Additional measurements were made on Novem- 
ber 23, December 1-2, and December 10, 2001. Our target 
list is given in Table Q Column 4 gives the nominal LSR 
velocity of each source adopted for the search. The source 
Ust consists of three well-studied regions of massive star for- 
mation w hich are known to be rich in complex organic mole- 
cules (e.g. lBlake et al.ll987HMiao et alll995l:lRemiian et alJ 
l2002l) . 

The heterodyne receivers used were the dual-channel 
B3 SIS mixers in the single-sideband mode (SSB) for the 
345 GHz submilliter band (315-373 GHz). The sideband 
rejection was ~13 dB (a factor of ^20). The SSB system 
temperatures in fair (medium) weather conditions - atmo- 
spheric opacity r(225 GHz) between 0.08 and 0.12 - were 
generally ~400 to 550 K. At 329.9 and 363.1 GHz, however, 
Tsys as high as ~700 K and ~900 K were recorded, due to 
poor atmospheric transmission and higher receiver tempera- 
ture. The half-power beamwidth (HPBW) of the telescope is 
~14 and the main-beam efficiency, T^mb, is 0.63. Data were 
obtained in the position-switching mode with offset 20 west 

^ The JCMT is operated by the Joint Astronomy Centre in Hilo, 
Hawaii on behalf of the present organizations: the Particle Physics 
and Astronomy Research Council in the United Kingdom, The 
National Raserch Coucnil of Canada and the Netherlands Orga- 
nization for Scientific Research 



in azimuth. Pointing and focus was checked regularly at a 2- 
hour interval. The resultant spectra are on the antenna tem- 
perature scale, T^, which has been corrected from chopper 
wheels calibration for atmospheric transmission and losses 
associated with rearward scattering. Since the spatial extent 
of each emission source is not known, no main beam correc- 
tion is applied. Further corrections for the forward scatter- 
ing and spillover efficiency (?7/ss = 0.82) convert the source 
antenna temperature, T^, to the source brightness temper- 
ature, T^. 

By employing the Dutch Autocorrelation Spectrometer 
(DAS) backend with two subsystems and a bandwidth of 500 
MHz for dual-polarity operation, we have a spectral resolu- 
tion of 756 kHz and a channel spacing of 625 kHz. Typically 
an integration time of 3 to 4 hours was achieved. Velocity 
shifting of ±4.5 km s~^ with respect to the nominal LSR 
velocity, about ±5 MHz with respect to the rest frequency 
in the 345 GHz band, was executed as a common practice 
during observations in order to neutralize the effect of pos- 
sible low-level gain variations in the DAS and to identify 
potential interlopers from the image sideband. The JCMT 
data were reduced using the Specx spectral line reduction 
package. 

Observations of submillimeter high-frequency transi- 
tions yield smaller telescope beams, preferentially sample 
the warmer and denser regions, and help to avoid line con- 
fusion with interloping emission from cooler envelope mate- 
rial along the line-of-sight. At hot core temperatures of 100 
K or more, searches in the submillimeter regime are further 
favoured by the fact that the higher frequency transitions 
are expected to be intrinsically stronger. 

Good agreement between observed line frequencies and 
those measured in the laboratory, for four or more spectral 
lines, are the minimum requirements for claiming identifica- 
tion of a new interstellar molecule. The rotational spectrum 
of pyrimidine has recently been measured over the spec- 
tral range 3-337 GHz, and the calculat ed dipole momen t 
of pyrimidine is /itot — = 2.39 Debye jKisiel et alJll999ll . 
Hence only fo-type transitions are observable. Our astronom- 
ical search was based on the best candidate transitions of 
pyrimidine in spectral regions free from known spectral line 
contamination. Four bandheads made up of closely spaced 
high-J transitions (at higher energy levels) plus 2 pairs of 
low-J doublet lines were observed in a total of 6 different 
spectral bands. The observed transitions are listed in Ta- 
ble 2; transitions with line strengths smaller than 10.0 are 
not included. Column 1 lists the line number; each line may 
include multiple pyrimidine transitions which are blended 
into one unresolved, single spectral-line feature. 

Figure Q illustrates the predicted relative line intensity 
in an arbitrary scale of the pyrimidine spectra over a fre- 
quency range of to 500 GHz at various rotational tem- 
peratures (Trot = 50 K, the top panel; 100 K, the mid- 
dle panel; and 200 K, the bottom panel). All possible spec- 
tral blends are accounted for in the simulation, which is 
made for a Gaussian lineshape with an assumed equivalent 
linewidth (FWHM) of 7 km s~\ Note that the linewidth, 
when Doppler converted from velocity to frequency, in- 
creases with frequency, which is conducive to the formation 
of stronger blends at higher frequencies. 



Table 2. Pyrimidine transitions observed. 



Line Rest Frequency Transition S*j Ej 



(MHz) 




- 








(cm" 




329961, 


.004 


53 1,53 


— 


52o,52 


52, 


.48 


288, 


,87 


329961, 


.004 


53o,53 


— 


52l,52 


52, 


.48 


288, 


,87 


329963, 


.042 


522,51 


— 


51l,50 


50, 


.45 


288, 


,46 


329963, 


.042 


52i,5i 


— 


512,50 


50, 


.45 


288, 


,46 


329963, 


.507 


43ll,33 


— 


42l0,32 


32 


.73 


266, 


,24 


329963, 


.507 


43 10, 33 


— 


42ll,32 


32, 


.73 


266, 


.24 


329963, 


.584 


42l2,31 


— 


41ll,30 


30 


.83 


261, 


.71 


329963, 


.584 


42ii,3i 


— 


41l2,30 


30, 


.83 


261. 


.71 


329963, 


.825 


44l0,35 




- 439,34 


34 


.65 


270, 


.36 


329963 


.825 


449,35 - 




43l0,34 


34 


.65 


270, 


,36 


329964 


.312 


41l2,29 


— 


40i3,28 


28, 


.95 


256, 


,77 


329964, 


.312 


41l3,29 


— 


40l2,28 


28, 


.95 


256, 


,77 


329964, 


.345 


459,37 


— 


448,36 


36, 


.58 


274, 


,06 


329964, 


.345 


458,37 


— 


449,36 


36, 


.58 


274, 


,06 


329964 


.448 


513,49 




502,48 


48, 


.43 


287, 


.64 


329964, 


.448 


512,49 




503,48 


48, 


.43 


287. 


.64 


329964, 


.911 


468,39 


- 


457,38 


38, 


.52 


277, 


.36 


329964 


.911 


467,39 


— 


458,38 


38, 


.52 


277, 


.36 


329965 


.303 


503,47 


— 


494,46 


46, 


.43 


286, 


.41 


329965, 


.303 


504,47 


— 


493,46 


46, 


.43 


286, 


.41 


329965, 


.396 


477,41 


— 


466,40 


40, 


.48 


280, 


,24 


329965, 


.396 


47(i,4i 


— 


467,40 


40, 


.48 


280, 


,24 


329965, 


.687 


494,45 


— 


485,44 


44, 


.43 


284, 


.76 


329965, 


.687 


495,46 


— 


484,44 


44, 


.43 


284, 


.76 


329965, 


.688 


486,43 


— 


475,42 


42, 


.45 


282. 


.71 


329965 


.688 


485,43 


— 


476,42 


42, 


.45 


282, 


.71 


329966 


.047 


40l4,27 


— 


39 13, 26 


27 


.08 


251, 


.42 


329966, 


.047 


40l3,27 


— 


39l4,26 


27 


.08 


251, 


.42 


336125, 


.535 


54 1,54 


— 


53o,53 


53, 


.48 


299, 


,88 


336125, 


.535 


54o,54 




53l,53 


53, 


.48 


299, 


.88 


336127, 


.005 


43l2,32 


— 


42ii,3i 


31, 


.81 


272, 


.72 


336127, 


.005 


43ll,32 


— 


42l2,31 


31 


.81 


272, 


.72 


336127 


.169 


44io,34 


— 


43ll,33 


33 


.72 


277, 


.25 


336127 


.169 


44ll,34 


— 


43l0,33 


33 


.72 


277, 


.25 


336127 


.422 


42l2,30 




41l3,29 


29 


.92 


267, 


.78 


336127 


.422 


42l3,30 


— 


41l2,29 


29, 


.92 


267, 


.78 


336127 


.563 


532,52 


— 


52i,5i 


51, 


.45 


299, 


.47 


336127, 


.563 


53l,52 


— 


522,51 


51, 


.45 


299, 


,47 


336127, 


.680 


45l(),36 




- 449,35 


35, 


.63 


281, 


,37 


336127, 


.680 


459,36 ■ 




44l0,35 


35, 


.63 


281, 


,37 


336128 


.356 


468,38 


— 


459,37 


37, 


.57 


285, 


,07 


336128, 


.356 


469,38 


— 


458,37 


37, 


.57 


285, 


,07 


336128, 


.740 


41l4,28 


— 


40l3,27 


28, 


.05 


262, 


.43 


336128, 


.740 


41l3,28 


— 


40l4,27 


28, 


.05 


262, 


.43 


336128, 


.946 


523,50 




512,49 


49, 


.43 


298, 


.65 


336128 


.946 


522,50 




513,49 


49, 


.43 


298, 


.65 


336129, 


.052 


478,40 


- 


467,39 


39, 


.51 


288, 


,36 


336129, 


.052 


477,40 


— 


468,39 


39, 


.51 


288, 


,36 


336129, 


.643 


486,42 




477,41 


41 


.47 


291, 


,24 


336129, 


.643 


487,42 




476,41 


41, 


.47 


291, 


,24 


336129, 


.762 


513,48 




504,47 


47, 


.43 


297, 


,42 


336129, 


.762 


514,48 




503,47 


47, 


.43 


297, 


.42 


336130, 


.023 


496,44 




485,43 


43, 


.45 


293, 


.71 


336130, 


.023 


495,44 




486,43 


43, 


.45 


293. 


.71 


336130 


.092 


505,46 




494,45 


45, 


.43 


295, 


.77 


336130, 


.092 


504,46 




495,45 


45, 


.43 


295, 


,77 


336131, 


.417 


40l5,26 




39l4,25 


26, 


.20 


256, 


.66 


336131 


.417 


40l4,26 




39l5,26 


26, 


.20 


256, 


.66 


338017, 


.612 


2727,1 




2626,0 


24, 


.97 


146, 


.26 


338018, 


.346 


2727,0 




2626,1 


24, 


.97 


146, 


.26 
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Table 2. (continued). 



sst Frequency 


Transition 


s: 


* 

ul 


Et 


(MHz) 


JKa,Kc 


~ Jk' K' 

a' c 




(cm" 




342289, 


.902 


55i,55 


- 54o,54 


54, 


.48 


311, 


,09 


342289, 


.902 


55a, 55 


- 54l,54 


54, 


.48 


311, 


,09 


342290 


,286 


44i2,33 


- 43ll,32 


32, 


.80 


283, 


,93 


342290 


.286 


44il,33 


— 43l2,32 


32, 


.80 


283, 


,93 


342290, 


.410 


43l2,31 


— 42l3,30 


30, 


.91 


278, 


,99 


342290, 


.410 


43l3,31 


— 42i2,30 


30. 


.91 


278, 


.99 


342290 


.681 


45l0,36 


- 44ll,34 


34. 


.70 


288, 


.46 


342290 


.681 


45ll,36 


- 44l0,34 


34. 


.70 


288, 


.46 


342291, 


.342 


42l3,29 


— 41 14,28 


29, 


.03 


273, 


.64 


342291, 


.342 


42i4,29 


- 41i3,28 


29, 


.03 


273, 


,64 


342291 


.379 


469,37 


- 45l0,3B 


36, 


.62 


292, 


,58 


342291, 


.379 


46l0,37 


- 459,36 


36, 


.62 


292, 


,58 


342291 


.921 


54l,53 


- 532,52 


52, 


.45 


310, 


,68 


342291 


.921 


542,53 


- 53l,52 


52, 


.45 


310, 


,68 


342292, 


.209 


478,39 


— 469,38 


38. 


.56 


296, 


.28 


342292, 


.209 


479,39 


— 468,38 


38. 


.56 


296. 


.28 


342293 


.032 


488,41 


- 477,40 


40. 


.51 


299, 


.58 


342293 


.032 


487,41 


- 478,40 


40. 


.51 


299, 


.58 


342293 


.280 


533,51 


- 522,50 


50, 


.43 


309, 


.86 


342293, 


.280 


532,51 


— 52-) tin 


50, 


.43 


309, 


.86 


342293, 


.488 


41l4,27 


- 40l5,26 


27, 


,17 


267, 


,87 


342293, 


.488 


41l5,27 


- 40l4,26 


27 


,17 


267, 


,87 


342293 


.729 


496,43 


- 487,42 


42, 


.47 


302, 


,46 


342293, 


.729 


497,43 


— 486,42 


42. 


.47 


302, 


.46 


342294, 


.057 


523,49 


— 514,48 


48. 


.43 


308, 


.63 


342294, 


.057 


524^49 


— 513,48 


48. 


.43 


308, 


.63 


342294, 


.195 


505,46 


- 496,44 


44. 


.44 


304, 


.93 


342294, 


.195 


506,45 


- 495,44 


44. 


.44 


304. 


.93 


342294 


.334 


515,47 


— 504,46 


46. 


.43 


306, 


.98 


342294 


,334 


514,47 


- 505,46 


46, 


.43 


306, 


,98 


342297, 


.447 


40l5,25 


— 39l6,24 


25. 


.33 


261, 


.70 


342297, 


.447 


40l6,26 


— 39l6,24 


25. 


.33 


261, 


.70 


348453, 


.267 


44i2,32 


— 43l3,3l 


31. 


.89 


290. 


.41 


348453, 


.267 


44l3,32 


— 43l2,31 


31. 


.89 


290, 


.41 


348453, 


.422 


45l2,34 


- 44ii,33 


33. 


.78 


295, 


.35 


348453, 


.422 


45ll,34 


- 44i2,33 


33. 


.78 


295, 


.35 


348453 


.838 


43l4,30 


- 42i3,29 


30, 


.01 


285, 


.06 


348453, 


.838 


43l3,30 


— 42l4,29 


30, 


.01 


285, 


,06 


348454, 


.040 


46ll,36 


- 45l0,35 


35, 


.69 


299, 


,88 


348454, 


.040 


46l0,36 


- 45ii,35 


35, 


.69 


299, 


,88 


348454, 


.101 


56i,56 


- 55o,55 


55 


.48 


322, 


,51 


348454, 


.101 


56o,56 


- 55l,55 


55 


.48 


322, 


,51 


348454, 


.919 


47l0,38 


- 469,37 


o7. 


.61 


304, 


.00 


348454, 


.919 


479,38 • 


- 46l0,37 


37. 


.61 


304. 


.00 


348455, 


.498 


42l4,28 


- 41l6,27 


28. 


.15 


279, 


.29 


348455, 


.498 


42l5,28 


— 41l4,27 


28, 


.15 


279, 


.29 


348455 


,899 


489,40 


- 478,39 


39 


.55 


307, 


,70 


348455 


,899 


488,40 


- 479,39 


39 


.55 


307, 


,70 


348456, 


.112 


552,54 


- 54l,53 


53 


.45 


322, 


,10 


348456, 


.112 


55l,54 


- 542,53 


53 


.45 


322, 


,10 


348456, 


.848 


498,42 


- 487,41 


41 


.50 


310, 


,99 


348456, 


.848 


497,42 


- 488,41 


41, 


.50 


310, 


.99 


348457, 


.447 


542,62 


— 533,51 


51, 


.43 


321, 


.28 


348457, 


.447 


543,52 


- 532,51 


51, 


.43 


321, 


.28 


348457, 


.650 


506,44 


- 497,43 


43, 


.47 


313, 


.87 


348457, 


.650 


507,44 


- 496,43 


43, 


.47 


313, 


.87 


348458, 


.186 


533,50 


— 524,49 


49, 


.42 


320, 


,04 


348458, 


.186 


534,50 


— 523,49 


49, 


.42 


320, 


,04 


348458, 


.202 


515,46 


— 506,45 


45, 


.44 


316, 


,34 


348458, 


.202 


516,46 


- 505,45 


45, 


.44 


316, 


,34 


348458, 


.410 


524,48 


— 516,47 


47, 


.43 


318, 


.40 


348458, 


.410 


525,48 


— 514,47 


47, 


.43 


318, 


.40 
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Table 2. (continued). 


Line 


Rest Frequency 
(MHz) 


Transition 
JKa,Kc - Jk^.k^ 




(cm ■"-) 




348458.771 
348458.771 


41l5,26 — 40l6,25 
41l6,26 — 40l5,25 


26.30 
26.30 


273.12 
273.12 


6 


363098.770 
363099.057 


2929,1 — 2828,0 
2929,0 — 2828,1 


26.97 
26.97 


169.23 
169.23 



* Line Strength. 

t Lower energy level. 

The strongest lines arise from blends of many transi- 
tions in high-J bandheads. These bandheads become more 
compact and have more lines at higher frequencies. On 
the other hand, following Boltzmann population distribu- 
tion, intensities of individual transitions forming the spectral 
blends, or bands, maximize at appreciably lower frequencies 
(henceforth the Boltzmann maxima) than the maxima ap- 
parent in the spectral profiles in Figure Q At Trot = 100 
K, for example, the strongest individual band line, i.e., the 
Boltzmann mciximum, though not obvious in the spectrum 
(the middle panel of Figure 0, is near 195 GHz. Never- 
theless, because at such frequencies the band compression 
effect becomes dominant and compensates for the intensity 
decrease of transitions at frequencies beyond the Boltzmann 
maximum, the overall band profile reaches an intensity max- 
imum at the much higher frequency of ~300 GHz. Similarly 
at Trot = 50 K and 200K, the Boltzmann maxima are at 
^^135 GHz and ~285 GHz, respectively, while at higher fre- 
quencies where band compression dominates, the apparent 
maximum intensity features for a linewidth of 7 km s^^ are 
predicted to be at ~180 GHz and ~335 GHz, accordingly. 
The 324 GHz spectral blend, which protrudes above the gen- 
eral profile and is seen in the middle and bottom panels of 
Figure arises from accidental addition of some non-band 
lines to the spectral blend. The prevailing spectral features 
at 100 K and 200 K are compact high-J R-type bandheads 
of the type observed in our search (see Table 2). For an ex- 
citation temperature near 100 K, it is clear that searches in 
the spectral region near 1-mm (300 GHz) should provide the 
best chance for a detection. 



3 RESULTS 

Table |21 lists the important measured physical parameters 
of all the pyrimidine lines observed. The pyrimidine lines 
are separated into three sources according to the actual ob- 
servations (Column 1). Column 2 gives the line numbers, 
which correspond to the line numbers shown in Column 1 of 
Table 2. Most of the pyrimidine lines are only observed in 
one or two sources but not in all. The fact that most of the 
data are fragmentary is mainly due to scheduling difficulties 
at the JCMT, leading to limited or uncertain availability of 
observing time for different sources. 

In total, 6 different pyrimidine lines were observed to- 
ward the 3 target sources. The peak antenna temperature of 
the spectral line observed, i.e. the T^ upper limit for a non- 
dection, is given in Column 3 in mK; "rms" indicates the l-u 
noise level per channel of the spectrum, where a = T\{rms). 



Column 4 lists the upper limits of the total column density 
Nf ot derived. T^, hence also Ntot, are not listed if the target 
pyrimidine line is completely overwhelmed by nearby strong 
interlopers. Self-explanatory comments for each spectral line 
observed in each source are given, when necessary, in Col- 
umn 5. Interlopers, which are largely blended with (denoted 
as "(B)") and smeared out the target pyrimidine lines or 
partially blended (denoted as "(PB)") with the candidate 
lines, are also included in Column 5 in the format mole- 
cule/frequency (MHz). 

Assuming that the pyrimidine lines are optically thin, in 
LTE, and that the rotational excitation temperature. Trot, 
is much higher than the background brightness temperature, 
the beam-averaged total column density of C-C4H4N2 can be 
written as: 



iVtot(cm-^) = (1.669x10^^)1^ 



Si,; fJ-l 



exp 



where 



Wk= Tr* dv = 



Ta* 



dv 



(2) 



in (K km s~^) is the integrated intensity of the spectral 
line either from a single pyrimidine transition, or from un- 
resolved multi-transitions such as in a bandhead. The rest 
frequency (MHz) is 1/, nt the permanent 6-dipole moment in 
Debye, Sui the line strength, E„ the upper energy level in 
K and Qrot the rotational partition function. All pyrimidine 
lines reported here contain more than one pyrimidine tran- 
sition, a summation of (SuiM5)/[exp(Eu/Trot)] over all rele- 
vant transitions was thus applied. In addition, a rotational 
temperature Trot = 100 K was assumed for all transitions. In 
the case of a non-detection, for an interloper-blended pyrim- 
idine line, is obtained from the T^ upper limit listed in 
Column 3 of Table |H] corrected for rifss- If the equivalent 
linewidth could not be unambiguously defined from the ob- 
served spectrum, the column density limits were evaluated 
using a representative value of the equivalent linewidth; a 
value of Av ~ 10 km s~^ was adopted for all three target 
sources. 

Sample pyrimidine spectra are shown in Figures |21 and 
13 With a limited number of spectral lines observed in each 
source, our submillimeter-wave search for pyrimidine did not 
yield a definite detection of pyrimidine in either Sgr B2(N) or 
Orion KL. By averaging over all observed lines measurable 
in each source (see Table , the inferred upper limits on 
the total column density in Sgr B2 and Orion are 1.7x10^* 
cm^'^ and 2.4x10^'* cm~^, respectively. 

In the case of W51 el/e2, there is a spectral feature at 
~ 57 km s~^ with respect to the rest frequency adopted for 
the observation (see the spectrum shown in the lower panel 
of Figure |3J. This feature coincides precisely with the J ^ 
55 band-head of pyrimidine between 342289.9 and 342297.4 
MHz (Line 4), containing 32 transitions at Vlsr = 59.0 km 
s^^. An unidentified U-line at 342290.0 MHz was previousl y 
reported in a line survey of Orion KL dSchilke. et al.|[r997l) : 
unfortunately we did not observe either Orion or Sgr B2 at 
this particular frequency. 

The total pyrimidine column density in W51 el/e2 de- 
rived solely from Line 4 is 4.8x10^^* cm~^; the averaged up- 
per limit to the column density from all 3 lines available in 
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Figure 1. The predicted relative line intensity of the pyrimidine spectra over a frequency range of to 500 GHz at various excitation 
temperatures (50 K, top; 100 K, middle; and 200 K, bottom). The ordinate is in an arbitrary scale. These spectra are a<;tually contour 
plots with point spacing of 0.5 MHz but have the appearance of stick diagrams because line contours are very narrow relative to the total 
breadth of the plotted spectrum. Each such stick is thus in fact a very compressed line profile for a spectral blend, and each spectrum 
is in principle 1000000 points long. The simulation is for a Gaussian lincshapc with an assumed equivalent linewidth (FWHM) of 7 km 
s-l. The apparent maximum intensity features are at ~180 GHz for 50 K, ~300 GHz for 100 K and ~335 GHz for 200 K. The 324 GHz 
spectral peak (middle & bottom) is due to accidental addition of some non-band lines to the spectral blend. 



W51 is Ntot ^ 3.4x10 cm" . It is interesting to note that 
the pyrimidine column density computed only from Line 4 
(assuming Trot = 100 K) is just slightly higher (within a fac- 
tor of 2) than the upper hmits deduced from the two other, 
interloper-contaminated, lines. It appears that if Line 4 of 
W51 is truly a detection, we may have expected to detect 
another bandhcad. Line 5, whose transitions are at iiigher 
frequencies. It could be argued that the fact that we failed 
to detect Line 5 makes the identification of Line 4 less likely. 



However, Line 5 is expected to be much weaker than Line 
4 at an excitation temperature lower than 100 K, and so, 

in the absence of accurately measured pyrimidine rotational 
temperatures for these sources, we cannot entirely rule out 
a true detection in W51. 

Finally, the submillimeter transitions targeted in this 
search possess rather high excitation levels, and so any ob- 
servable pyrimidine emission would mostly originate from 
regions of high temperature and/or density, if radiative ex- 
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Table 3. The measured physical parameters of the pyrimidine lines observed. 



Source 


Line 


± rms 


Ntot 


Comment 






(mK) 


(10" cm-2) 




Sgr B2(N) 


1 






(B) C2H5OH/329954.9; (B) C2H5OH/329956.5; (B) C2H5OH/329968.2 


Sgr B2(N) 


3 


52.6±13.5 


<1.68 


(B) CH3OCH3/338O25.5 


Sgr B2(N) 


5 






(B) CH3CHO/348449.6; (B) C2H5CN/348473 


Orion KL 


2 


215.7±13.8 


<1.18 


(B) C-C3H2/336I28.5; (B) C2H3CN/336I37.8; 










(PB) HCOOCH3/336III.3; (PB) SO2/336II3.5 


Orion KL 


5 


400.0±12.6 


<3.66 


contaminated by a spectral feature from the image sideband 


Orion KL 


6 






(B) CH3C3N/363097.I; (B) Sii8o/363100.7; (B) C2H5CN/363IO7 


W51 el/e2 


3 


140.6±12.9 


<3.14 


(B) CH3OCH3/338O25.5 


W51 el/e2 


4 


631.3±23.6 


<4.79 


detected? (U-342290); (PB) C2H3CN/342286.8 


W51 el/e2 


5 


127.5±12.3 


<2.20 


(B) CH3CHO/348449.6 



citation is negligible. As a result, these excited pyrimidine 
molecules would have been located in regions fairly close to 
protostellar hot cores within molecular clouds. It is therefore 
probable that the real source size is smaller than the tele- 
scope beam (^^14 ); it is more likely to be the case particu- 
larly for distant sources such as Sgr B2(N) and W51 el/e2. 
Hence, one should note that our observations would suffer 
from beam dilution if the source extent is indeed smaller 
than the beam. The derived upper limits on the column 
densities, based on the assumption that the source fills the 
beam, could consequently be underestimated. 

The upper limits of fractional abundance of pyrimi- 
dine with respect to molecular hydrogen, X(c— C4H4N2) 
~ '^tot{c— C4H4N2)/Ntot(H2), may also be deduced. Molec- 
ular hydrogen column densities inferred from single-dish ob- 
servations with beam sizes similar to the JCMT are em- 
ployed, in order to determine the beam-averaged abundances 
more accurately. The H2 column densities adopted arc: ~ 5 
X 10^" cm~^ for Sgr B2(N) (Numniclin ct al. 2000), 8 x 
10^^ cm"^ for Orion K L llSutton et a l. 1995), and 1 x 10^ ^ 
cm"^ for W51 el/e2 ij.Taffe. Becklin fc Hildebrandl llQSA . 
The upper limits of pyrimidine fractional abundances thus 
estimated are X(c - C4H4N2) ^ 3.4 x 10"" for Sgr B2, 3.0 
X 10"^° for Orion, and 3.4 x 10"^" for W51. 



4 DISCUSSION 

Our submillimeter search for the nucleic acid building-block 
pyrimidine has not been successful. We have perhaps de- 
tected one-/jne in a single source but this needs to be con- 
firmed. The negative result therefore cannot be considered 
definitive and searches at lower frequencies may yet detect 
pyrimidine in hot cores. The low inferred pyrimidine abun- 
dance limits may simply reflect the fact that pyrimidine is 
generally of low abundance in massive star-forming cores, 
and hence it is difficult to populate the high-J energy lev- 
els (J ^ 56 in our study) which are observable at submil- 
limeter wavelengths. On the other hand, it could also imply 
that either pyrimidine is not evaporated efficiently from dust 
grains at hot core temperatures, or is destroyed easily in the 
warmest regions in hot cores sampled by submillimeter ob- 
servations. 

It would appear that the best chance for an astronom- 



ical detection of pyrimidine, and other Nitrogen Heterocy- 
cles, is probably in the circumstellar envelopes (CSEs) of 
carbon stars. AGB and post-AGB stars (e.g. IRC-l-10216, 
CRL 618 and CRL 2688) are copious producers of carbona- 
ceous dust particles. Recent ISO observations of well-known 
protoplanetary nebulae (PPNe) have uncovered more new 
organic molec ules, including the first detectio n of benzene 
in CRL 618 JCernicharo et ahlbOOlal l2001b^ . The initial 
stages of dust formation involve polymerization of acety- 
lene to form benzene and subsequent C2H2 additions lead 
to la rge polycyclic aromatic hydrocarbo n (PAH) molecules 
(e.g. ICherchneff. Barker fc Tielen!Jll^92^ . During this reac- 
tion sequence other triply-bonded molecules can also add 
to the growing ring structures. In particular, recent theo- 
retical work indicates that N atoms can become incorpo- 
rated in ring structures through a dditions involving HCN 
jRicci. Bauschlicher fc Bakes! l200lD . However, the kinetics 
of ring growth suggests that single rings containing two N 
atoms will be less favoured than single rings with one N 
atom (i.e. pyridine). Even two-ring compounds containing 
a single N atom, such as quinoline and isoquinoline, the N- 
substituted analogues of naphthalene, co uld be more abun- 
dant than pyrimidine jKisiel et alJl2003l) and may explain 
why our search appears to have been unsuccessful. Searches 
for pyridine, quinoline and isoquinoline in the molecular en- 
velopes of evolved stars are currently underway. 

As an alternative to r ing fo rmation in neutral-neutral 
reactions, IWoods. et al.l ll2002f) have shown that ion- 
molecule reactions are able to account for the abundance 
of benzene observed in CRL 618. The growth of larger 
rings, and incorporation of heteroatoms into them, has not 
yet been consi dered in circumstellar ion-molecule chemistry. 
IWoods. et alJ (|^003) have proposed one neutral process, in- 
volving benzene and CN, that could lead to nitrogen being 
added to ring structures. However, the product of this reac- 
tion has the nitrogen present in a side-group (i.e. benzoni- 
trile, c-CeHsCN) and not bonded into the ring. 

Another possible reason for a low interstellar pyrimidine 
abundance concerns its photostability. The infrared spectra 
and photostability of pyrimidine have recently been mea- 
sured in an Ar matrix at 10 K dPeeter et al J 12003). The 
stability of pyrimidine (which does not have real aromatic 
properties) against UV photolysis is rather limited and far 
below that of mono and polycyclic aromatic rings. These ex- 
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Figure 2. Sample pyrimidine spectra at 338.0 GHz (Line 3; 
2727,1-2626,0 and 2727,o-2626,l) of Sgr B2(N) (top) and of W51 
el/e2 (bottom). The two vertical lines mark the two expected 
transitions of pyrimidine. The abscissas give the LSR velocities 
with respect to the rest frequencies adopted for the observations 
at the nominal LSR velocities 64.0 km for Sgr B2(N) and 
60.0 km s-i for W51 el/e2. 

periments appear to indicate that any circumstellar pyrim- 
idine would be easily destroyed by UV photons, or cosmic 
ray particles, soon after delivery to the interstellar medium. 
In this case, any interstellar pyrimidine would have to be 
produced in dark molecular clouds. 



5 CONCLUSION 

To strengthen the role of prebiotic interstellar matter in As- 
trobiology, we have searched for 6 pyrimidine lines in three 
massive star-forming regions: Sgr B2(N), Orion KL and W51 
el/e2. Our search was unsuccessful and did not yield a con- 
clusive result, with only one potential single-Zine detection. 
The abundance limits inferred are ^ 3.4 X 10"" for Sgr 
B2, 3.0 X 10"^° for Orion, and 3.4 x 10"^° for W51. Our 
negative result may simply reflect the fact that interstellar 
pyrimidine is of low abundance. Sources that are ongoing 
sites of carbon dust formation, such as the C-rich envelopes 
of AGB and post- AGE stars, probably present the best op- 
portunity for detecting pyrimidine and other nitrogen hete- 
rocycles. 
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